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Abstract 

We study the influence of penguin (especially, electroweak penguin) effects 
on some methods of measuring the angles a, j3, and 7 in the CKM unitarity 
triangle. We use next-to-leading order effective Hamiltonian, and present 
numerical estimates based on the factorization approximation. We find that 



in 
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some techniques suggested in the literature, especially for a determination, 

are not workable in light of the electroweak penguin effects. Nevertheless, 

there are methods that would work for each angle determination. For angle 

(3 we consider B — > D + D~ mode and estimate the penguin contamination. 

For angle 7 we consider a method based on SU(3) symmetry and carefully 

consider SU(3) breaking effects. We point out regions in the parameter space 

where this method could be used reliably. 
PACS numbers:11.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 



The measurement of the e-parameter in the K° — K° meson system is the only direct 
evidence for CP violation in the laboratory Many models have been proposed to explain 
this phenomena The Standard Model (SM) of three generations with the source for 

CP violation arising from the phases in the Cabibbo-Kobayashi-Maskawa (CKM) matrix is 
consistent with the experiment M. It is necessary to perform more experiments to find out 
the source or sources of CP violation and to test the CKM model. A unique feature of the 
CKM model for CP violation is that the CKM matrix is a 3 x 3 unitary matrix. Due to 
the unitarity property, when summed over the row or column of matrix elements V!y times 
complex conjugate matrix elements V* k , the following equation holds, 

i 

This equation defines a triangle when j ^ k. For example, for j = d and k = b a triangle 
shown in Fig. 1 with three angles a = Arg(-V td V t * b /V* b V ud ), (3 = ATg(-V cd V^/V t tV td ), and 
7 = Kig{—V ud V* b /V* b V cd ) is defined. These angles are related to phases of the CKM matrix 
elements. In Wolfenstein parametrization the three angles of the triangle are given by a = 
Arg (-V td /V* b ), (3 = Arg V t * d , 7 = Arg V* b 0. The sum of these three angles must be 
equal to 180°. This is a unique feature of the CKM model for CP violation. This property 
provides an important way to check the validity of the CKM model if enough independent 
measurements of the sides and angles of the triangle can be performed experimentally. B 
meson decays provide a fertile ground to carry out such a test ||. 

Many methods have been suggested for measuring a, [3 and 7 using B decays |5|-[T0 



One class of methods involve the measurements of CP asymmetries in time evolution of B° 
decays into CP eigenstates. Such methods make it possible to measure the three angles of the 
unitarity triangle independently and without hadronic uncertainties if amplitudes depending 
on a single CKM phase dominate the decay process. For instance 0, sin 2a, sin 2(3, and 
sin 27 can be measured in decays B° — > 7r + 7r~, B° — ► ^K$, and B® — > pK$, respectively. 
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However, in practice due to simultaneous contributions to these decays from the tree and 
loop (penguin) effects, the measurements are much more difficult. In many cases the CKM 
phases can not be extracted just from these asymmetries. Additional information is needed. 
Some methods have been developed to extract CKM phases using relations based on isospin 
or flavor SU(3) symmetries to isolate the CKM phase of one type of amplitude such that 
the CKM phase determined is not contaminated by the presence of other amplitudes ||7Hl6|. 
Some of these relations when used for charged B meson decay modes, also provide new 
methods to measure some of the CKM phase angles. 

There are two types of penguin contributions, the strong and the electroweak penguins. 
Naively, one would expect that the electroweak penguin effects are suppressed by a factor of 
otem/ois compared with the strong penguin effects, and therefore can be negelected. Many 
previous methods for measuring the CKM phases have explicitly made such assumption. In 



a recent paper by two of us |12j , we pointed out that this assumption turns out to be wrong 
for large top quark mass. Some of the methods proposed in the literature become invalid 
when the electroweak penguin effects are included. In this paper we study the effects of the 
electroweak penguin on several methods for measuring the CKM phases in the literature. 
We will concentrate on methods based on B decays. There are methods based on B s decays, 
which are much more difficult to perform experimentally, and we will not discuss them in 
this paper. 

The paper is organized as following: In section II we present the full effective Hamiltonian 
responsible for B decays, and some isospin and SU(3) analysis of the decay amplitudes; In 
section III, IV and V we study the influence of penguin effects on some methods for measuring 
the CKM phases a, (3 and 7, respectively; And in section VI we present our conclusions. 

II. THE EFFECTIVE HAMILTONIAN FOR B MESON DECAYS 

The effective Hamiltonian up to one loop level in electroweak interaction for hadronic B 
decays can be written as 
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Hab=i = -7=-[V ub V^( Cl O^ + c 2 O u 2 ) + V A V^c x O\ + c 2 2 c ) - V^V* £ cA] + F.C, (2) 

V ^ i=3 

where Oj's are defined as 

0{ = q^Lfpfp-fLba, 0{ = gj^Lff-fLb, 

3( 5) = qi^Lb^rL(R)q\ O m = q alfl LbpZqp^L(R)q' a , 

3 3 
07(9) = -^qi^Lbi:e ql q'YR{L)q\ 08(io) = -q a -f^ 1 LbpEe q/ q'pj ,x R(L)q' a , 

On = JjL-m b q<j» v RT a bG>i v , Q 12 = -^-m^RbF^ , (3) 

where = (1 =F 7s)/2, / can be u or c quark, g can be d or s quark, and g' is summed 

over u, d, s, and c quarks, a and (3 are the color indices. T a is the SU(3) generator 
with the normalization Tr(T a T b ) = S ab /2. G^ u and F^ v are the gluon and photon field 
strength, respectively. c$ are the Wilson Coefficients (WC). Oi, 2 are the tree level and 
QCD corrected operators. O3-Q are the gluon induced strong penguin operators. are 
the electroweak penguin operators due to 7 and Z exchange, and "box" diagrams at loop 
level. The operators Oii,i2 are the dipole penguin operators. 

The WC's q at a particular scale /1 are obtained by first calculating the WC's at mw 
scale and then using the renormalization group equation to evolve them to fi. We have 
carryed out this analysis using the next-to-leading order QCD corrected WC's following 
Ref. Hl7j . Using a s {mz) = 0.118, a em {mz) = 1/128, m t = 176 GeV and /x « m& = 5 GeV, 
we obtain from top-quark contribution |l8j 

ci = -0.3125, c 2 = 1.1502, c 3 = 0.0174, c 4 = -0.0373, 
c 5 = 0.0104, c 6 = -0.0459, c 7 = -1.050 x 10~ 5 , 

c 8 = 3.839 x 10~ 4 , c 9 = -0.0101, ci = 1.959 x 10~ 3 . (4) 

It is interesting to note that the coefficient eg arising from electroweak penguin comtribution 
is not much smaller than coefficients of the strong penguin. This enhancement is caused 
by a term in the electroweak penguin contributions in which the WC is proportional to the 
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square of the top quark mass due to Z exchange. In some application we will need absorptive 
parts of c and u loop contributions. These are given in Ref. ||18|| . 

The coefficients for the dipole penguin operators at the two loop level have the following 



values fl9l| : 



cn = -0.299 , C12 = -0.634 . (5) 

To study exclusive B decays, we also need to transform the quark operators into hadrons. 
This is a difficult task. At present there is no reliable way to carry out this calculation. 
Nevertheless, many models and suggestions have been made to provide some handle on 
the related hadronic matrix elements. Symmetry considerations provide very powerful con- 
straints on the matrix elements and relate different decay amplitudes. Isospin and flavor 
SU(3) symmetry are two very useful symmetries used in the analysis in this paper. 

We can always parametrize the decay amplitude of B that arises from quark subprocess 
b — > uuq as 



A =< final state\Hl ff \B >= V uh V^T{q) + V tb V*P(q) , 



(6) 



where T(q) contains the tree as well as penguin contributions, while P(q) contains purely 
penguin contributions. 

When q is fixed, isospin symmetry relates some of the decay amplitudes generated by 
the effective Hamiltonian. In the case for q = d, isospin symmetry relates decay amplitudes 
for different B — > 7T7T or B — > pn decays. It also gives information on which operators in the 
effective Homiltonian contribute to certain isospin decay amplitudes. This is very important 
for our discussions in the rest of the paper. 

For q = d, the tree operators Oi i2 and the electroweak penguin operators O 7 _i contain 
AI = 1/2 and AI = 3/2 interactions whereas the strong operators Oz_§ and the dipole 
penguin operators On ; i2 contain only AI = 1/2 interaction. In the case of B — > titt, Bose 
symmetry requires nir to be in / = or / = 2 state. Since B meson is a /= 1/2 state, we 
immediately know that the strong and dipole penguin operators will not contribute to / = 2 



decay amplitude, for example B~ — > 7r~7r°. In the case for B — > pre, the final states can 
have / = 0, 1 or 2 states. We also easily see that the strong and dipole penguin operators 
will not contribute to the 1 = 2 decay amplitude. 

For q = s. The tree operators and the electroweak penguin operators Oj-iq contain 
AI = and 1, and the strong and dipole penguin operators O^q and On^2 contain AI = 
interaction only. In B — > ixK decays, the combinations, A{B~ — > ti K~) + A(B° — > 7i°K°) 
and A(B° — > tt + K~) — A(B~ — > tt~K°) contain only / = 3/2. We immediately know that 
the strong and dipole penguin will not contribute to these amplitudes. 

Isospin symmetry will not relate the amplitudes between the amplitudes with q = d 
and amplitudes with q = s. However, if one enlarges the symmetry group to the flavor 
SU(3), these can be related. The isospin relations will still be maintaned because isospin is 
a subgroup of flavor SU(3) symmetry. We shall now use SU(3) symmetry to obtain some 
relations which will be used in our later discussions. 

SU(3) relations for B decays have been studied by several authors |T3| , p3| , p0| , pT| . The 



operators Oi )2 , 3 _ 6> ii ; i2, and O 7 _i transform under SU(3) symmetry as 3 a + 3b + 6 + 15, 
3, and 3 a + 3;, + 6 + 15, respectively. In general, we can write the SU(3) invariant amplitude 
for B decay to two octet pseudoscalar mesons. For the T amplitude, for example, we have 

T = Af^B i H(3) i (Mj"M l k ) + C^A^^y 
+ Al^H^fi MjMf + Cl^M^QtMl 

+ Aj—^BiH (15 ) jf M l j Mi + B i M i i H(JF)\ k M l k , (7) 

where Bi = (B~,B°,B®) is a SU(3) triplet, M\ is the SU(3) pseudoscalar octet, and the 
matrices H represent the transformation properties of the operators 0\-\2- H(Q) is a traceless 
tensor that is antisymmetric on its upper indices, and H(15) is also a traceless tensor but 
is symmetric on its upper indices. We can easily see that the strong and dipole penguin 
operators only contribute to A 3 and C3. 

For q = d, the non-zero entries of the H matrices are given by 

iJ(3) 2 = l, H(6)l 2 = H(6)f = 1 , H(6)f = H(6)f = -1 , 
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tf(15)} 2 = #(15) 21 = 3 , H{lb)f = -2 , H(15)f = if (15)f = -1 . (8) 

For q = s, the non-zero entries are 

H(3f = l, H(Q)\ 3 = H(Q)f = l, H(6)f = H(Q)f = -l, 
H(T5)? = H(T5)f = 3 , H(15)f = -2 , H(15)f = HQS)? = ~ l ■ (9) 

In terms of the SU(3) invariant amplitudes, the decay amplitudes T(nn), T(irK) for B° — > 
7T7T, B° — > nK are given by 

T(n+n-) = 2Af 5) + C ( T 3) - Af 6) + Cf 6) + + 3C ( ^ , 
T(7r°7r°) = -^=(2A(3) + _ 4) + ^) + ^(15) ~~ 5( ^(T5)) ' 
T ("°) = ^7f ^(Ts) ' 

T(n-R°) = Cg } + Af 6) - Cg } + 3A^ } - , 

= + AJ 6) - Cf 6) + 3^ + 7C^) , 

T(vr + K-) = Cg) - Af 6) + Cf 6) - A T m + SCf^ , 



T(vr K ) - - ^/=(C(3) _ 4) + C(6) _ ^(15) _ 5( %5)) ' 
= ^(-C® - 4) + " 34l5) + 9C, (T5)) > 



We also have similar relations for the amplitude P(q). The corresponding SU(3) invariant 
amplitudes will be denoted by Af and Cf. It is easy to obtain the following relations from 
above: 

V2A(B° -> 7T 7T ) + V2A(B~ -> 7T-7T ) = 1(5° -> TT+O , 

v^A(S" -> - 2A(S" -> 7T-i?°) = >/fU(£~ -> r] 8 K-) . (11) 

One expects the hadronic matrix elements arising from quark operators to be the same 
order of magnitudes, the relative strength of the amplitudes T and P are predominantly 
determined by their corresponding WC's in the effective Hamiltonian. However, in order 



to numerically compare contributions from different operators, we have to rely on model 
calculations. In our later analysis, when such numerical calculations are required, we will 
use factorization approximation. These numerical numbers may not be accurate, but they 
will serve well in providing an idea of the validity of certain assumptions made. 

III. MEASUREMENT OF THE PHASE ANGLE a 

In this section we study the electroweak penguin effects on several methods proposed for 
measuring the CKM phase angle a. 

(1) From time dependent asymmetries in B — > tttt decays 

Let us first consider the standard method to measure the CKM phase a in B —>■ titt 
IK?]]. We present it here to set up our notations and also to clarify some issues. The time- 
dependent rate for initially pure B° or B° states to decay into a final CP eigenstate, for 
example 7r + 7r~ at time t is H 



r(B°(t) -> 7T+7T-) = \A\ 2 e- n [ - + ^ + 1 J A| cos (AMt) - ImX sin (AMt)], 
T(B°(t) -> 7T+7T-) = |A|V rt [ 1+ J A|2 - 1 ~J A|2 cos (AMt) + ImX sin (AMt)], (12) 



2 



where A is defined as 

A ^ 9 4 (13) 
with A = A(B° — > 7T + 7r _ ) and A = A(B° — > 7t + 7t~). Here p and q are given by the relations 

\B L , H >= p\B° > ±q\B° > (14) 
for the two mass eigenstates Bh and Bl. In the SM, 

g _ v ttVtd (15) 
P v tb v t */ 1 ; 

The parameter ImX can be determined by mesuring the coefficient in CP asymmetry in 
time evolution varying with time as a sine function. 
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Using the effective Hamiltonian given in section II, we can parametrize the decay ampli- 
tude in general as 

A = V ub V: d T + V tb V t * d P . (16) 
If the penguin amplitude can be neglected, we have 

ImX = Im(Yf^^Yk) = sin(2 a) . (17) 

y v tb v td v: b vj y } K } 

The phase a can therefore be determined. Similarly one can obtain a from B° — > 7r°7r°. 
When the pegnuin effects are included, ImX is not equal to sin(2a) any more. One finds 



2cos5sina + sin(2a)(i?-2cos(5 + a)) 

Asm (2a) = ImX — sin (2a) = —R — , 18 

V ' K 1 1 + R 2 - 2Rcos(5 + a) ' V ; 

where R = \P/T\, and 5 is the relative strong rescattering phase between the T and P 



amplitudes. It was estimated in Ref. [12] using factorization approximation that R = 7%. 
Using this number it was found that the error on the determination of the phase a can be 
as large as 12°. The error is even larger if B° — > 7r°7r° is used, where R is estimated to be 
23%. 

(2) Gronau-London method 

When penguin effects are included, the parameter ImX can be parametrized as 

ImX = |4. S i n (2 a + Q) . (19) 

The ratio can be determined by measuring the coefficient of the term in CP asym- 

metry in time evolution varying as cosine function at asymmetric colliders, and also at 
symmetric colliders p2" |. If 9 can be determined independently, the phase a can also be 



determined. To determine 6, Gronau and London M proposed using isospin relation 

V2A(B° -> ttV) + V2A(B- -> tt-tt ) = A(B° -> tt+tt-) , (20) 

obtained in eq. (|TTJ) , and similar relation for the CP-conjugate amplitudes for the corre- 
sponding anti-particle decays. If all the six amplitudes can be measured, the angle 6 can 
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be determined up to two fold ambiguity as shown in Fig.2. In this method GL argued, 
because B~ — > tt^tt has only 1 = 2 amplitude in the final state to which only AI = 3/2 
interaction in the Hamiltonian contribute, that A(B~ — ► 7r~7r°) and A(B + — > 7r + 7r ) have 
no contribution from penguin operators and therefore \A(B~ — > 7t~tc°)\ = \A(B + — > 7r + 7r°)| 
. The strong penguin only has A J = 1/2 interaction, so the strong penguin does not 
contribute to this decay. If the electroweak penguin effects are neglected, the equlity 
\A(B~ — > 7r~7r°)| = \A(B + — > 7r + 7r°)| is exact. The electroweak penguin actually con- 
tains AI = 3/2 interaction, and contributes to the decay amplitude. However, since the 
electroweak penguin < strong penguin < tree contribution for this process, the contribution 



is expected to be very small. An estimate based on factorization gives less than 3% |]T2 
The inclusion of the electroweak penguin effects can be safely neglected. This method, in 
principle, can determine the phase a at a few percent level. 
(3) Hamzaoui-Xing method 

A method to measure the phase a without using CP asymmetry in time evolution has 
also been proposed recently by Hamzaoui and Xing [|T1|] . We show that this method actually 



fails. Based on isospin consideration and factorization approximation, they parametrized 
the decay amplitudes for B — > tttt as follows: 

A +0 = A{B + - vr+vr ) = ~^Te« 

A + _ = A{B° -> 7T+7T-) = -Te i7 - P + _e i{& -P\ 

Ax. = A(B° -> vrV) = -^f<* + ^Pooe^\ (21) 

where the T and P^s are the tree and penguin amplitudes, respectively. 5 denotes the 
strong relative phase between the penguin and tree amplitudes. The parameter a denotes 
the color-mismatched suppressed contribution in the tree amplitudes. If electroweak penguin 
effects are neglected, 

P + _ = Poo = P- (22) 

One then obtains 
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cos(a + 5) = — [1 + x 2 ~ (1 + a) 2 R + _], 

cos(a - 5) = + x 2 - (1 + a) 2 i? + _], (23) 

where \ — P/T, and and _R + _ are the measurable quantities defined as R+j = 

\Aij\ 2 / (\A +0 \ 2 + |y4_ | 2 ) {A denotes the CP-conjugate amplitude of A). The parameters 
a and % can also be expressed in terms of experimental measurables, 

Roo — Roo 



a = -2 



X = yJ{l + a){aR+- + 2Roo)-a. (24) 

Therefore, if the assumptions are correct, the phase a could be determined. 

In order to obtain the above equations, a cruical assumption has been made that the 
parameter x+- — P+-/T is equal to Xoo — Poo/T. This equality is true only if there is no 
electroweak penguin contribution. The validity of the proposed method can be checked when 
electroweak penguin effects are included. Using factorization approximation, we obtain 

G F 

T = -—^\V ub V u d\{^ci + c 2 )T nn , 

Q 

P+- = -^\V ub V ud \ [£c 3 + c 4 + 2(£c 5 + c 6 + £c 7 + c 8 )X 1 + £c 9 + c X q]T^ , 

G 113 

Poo = -^|K&Kd| [£c 3 + c 4 + 2(£c 5 + c 6 -£c 7 - -c 8 )X 2 + -(c 7 + £c 8 - c 9 - £ci ) 

--(fo + cio)]^, (25) 
where £ = 1/JV with iV being the number of color, and 



(m b - m u )(m d + m u ) ' 2m d (m h - m d ) ' 

= ifAft«)(ml ~ ml) + fMX\ ■ (26) 



In our calculations we will use U = 132 MeV, and the form factors j% in Ref. @. For b 
quark mass, we use = 5 GeV, and for the light quark masses m u = 5.6 MeV, m d = 8.7 
MeV. We also treat £ as a free parameter and use the experimentally favoured value £ ~ 1/2 



24 1 . We find that for B° — > Tr°Tr° the ratio of the electroweak penguin to the strong penguin 
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amplitude is quite large (about 33%), while that ratio for B° — > ty + ti^ is only 5.7%. In 
B° — ► 7T°7r decay, the electroweak penguin contributions have the opposite sign to the 
strong penguin contributions and reduce the total penguin effects. We find the values for 
X+- and xoo are very different, x+-/Xoo — 1-71. Here we have neglected small contributions 
from the dipole penguin operators which do not affect the result significantly. It is clear that 



the assumption \+- — Xoo is badly violated, and therefore the method proposed in Ref. JTT 

fails completely. 

(4) Using B — > pn decays 

In Ref. ||,[| it was pointed out that neutral B — > pn can be used to extract a without 
ambiguities due to penguin contributions. It was shown in Ref. |8|U that by studying full 
Dalitz plot and time dependence for B° — > 7r + 7r~7r°, the amplitudes and phases of B° — > pn 
decays, S 3 = A{B° -> P + ^~), S± = A(B° -> p-vr+), 5 5 = -2A(5° -> A°), and their 
CP-conjugate amplitudes can be determined. Isospin analysis then shows that the sum 
S = S 3 + 5*4 + 5*5 has only 1 = 2 amplitude. Therefore it arises from AI = 3/2 interaction. 
If electroweak penguin effects are neglected, S has only tree contribution which contains 
the phase angle 7. Combined with angle j3 from the mixing parameter q/p in B°-B°, and 
considering coefficient of sin (AMt) (see eq.(|T^)), the phase a can be determined. 

Since the electroweak penguin contains AI = 3/2 interation, when its effects are included, 
S contains in addition to the tree amplitude S tree , also the electroweak penguin amplitude 
S ew which has a different weak phase. The determination of the phase a will therefore be 
contaminated. In the factorization approximation, we obtain 

s tree = + 0(ci + c 2 )(c + f ), 

S ew = ^V t lV td [3(^7 + c s )Xf - ^(c 7 + £c s )(C - f) 

-^(l+O(c 9 + c 10 )(C7 + T)], (27) 



where 



X 



2 

mi 



2m d (m b + m d ) 
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C = 2m p (e;.p^f p f+ 7r (m 2 p ), 

f = -2m p (e*-p 7T )f 7T A* p (ml). (28) 



Here we have neglected a small contribution due to annihilation effects. Note that the strong 
and dipole penguin operators do not contribute to S. For our numerical calculations we will 
use f p = 221 MeV, and the form factors f£ n and Aq P calculated in Ref. p3] . We find 

^4 « l-mV td \/\V ub \). (29) 

We see that the ratio of the electroweak penguin to the tree contribution is very small and 
S is dominated by the tree contribution. Therefore, this method of measuring the phase a 
is good to a few percent level. 
(5) Using B — > nK decays 

A method for measuring the phase angle a has also been proposed using |A>S| = 1 B 
decay processes by Nir and Quinn |J. Measurement of CP asymmetry in time evolution in 
B° — > n°Ks will be able to determine the parameter 

JmA = /m W(ffl^^ S ) )=/m(e A^^) ( 30 ) 

If penguin effects are neglected, A/A = e _2n and so JmA = sin(2a). For this decay it is 
obviously wrong to neglect the penguin effects because the penguin contribuitons are en- 



hanced by a factor of \V t bV t * s /V u bV* s \ ~ 50 |26| compared to the tree contributions. Even 



though the WC's of the penguin operators are smaller than the tree WC's, the net pen- 
guin contributions may be larger than the tree contributions. In Ref. || a method have 
been proposed to overcome difficulties associated with strong penguin effects by determin- 
ing e 2rr A(5° — > 7i°Ks)/A(B° — > ti Ks) directly from isospin analysis. This method requires 
measurement of the decay amplitudes, A(B° — > ir°K°), A(B° — >• 7r + K~), A(B~ -> n°K~), 
A(B~ — > 7i~K°), and their CP-conjugatey amplitudes. The strong and dipole penguin oper- 
ators do not contribute to the follwoing combinations because they are 1 = 3/2 amplitudes, 
which can also be easily seen from eq. ([II]) , 
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U EE ±A{B- -> 7T°K-) + ^(5° -> 7T°i? ) , 
P EE -> 7T-K°) + ±A(B° -> TT+i^-) , 

[/ = 1^(5+ -> 7T°K + ) + ^4(5° -> 7T°ir ) , 

V EE -> 7T + ^°) + ^(5° -> . (31) 

If the electroweak penguin effects are neglected, 

U = Ue i2 ^ , V = Ve i2 ~< . (32) 

When the eight decay amplitudes for B — > nK are measured, using the conditions in eq. , 
the quadrilaterals in eq. (O) 

A(B~ -> 7T ^-) - -^1(5- -> 7T-/?°) = 1(5° -> 7T /? ) + -^i(5° -> n + K-) , 
V 2 v 2 

A(B + -> 7T°K + ) - ^A(fi+ -> 7T + /^°) = A(5° -> 7T K°) + ^A(5° -> 7T-^+) , (33) 

v 2 y2 
can be constructed with U + V being a common diagonal. Here A's are defined as 
/!(£> — > 7riT) ee A(S — > 7riT)e 217 . Once these quadrilaterals are constructed, the quantity 
e l27 A(5° — > tt°Ks)/A(B° — > tt°Ks) can be easily determined. Combining this information 
with eq.(^), the phase a can be determined. 

However, since the electroweak penguin operators also contain AI = 1 interaction, we 
have to check the validity of the relations of eq.(^). In factorization approximation we find 
the magnitudes of the amplitudes U, V, U, and V: 

U = { OU2( Cl + ic 2 )C + (fo + c 2 )(T> + A')] 

- V t * b V ts [-3(c 7 + £c 8 - C9 - £c 10 )C' + 3(£c 7 + c 8 )(YT' + ZA') 
+ l^c 9 + c w )(T' + A)}}, 

v = %?t/2 { ^"Kfr + C2)(T ' " A) 

- V^V^fo + c 8 )(YT' - ZA) + ^(£c 9 + c 10 )(T' - A)]}, (34) 

where 
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(m b + m u )(m s - m u ) ' 
C = *fAfBK«)« - m 2 K ) + fsAOO, 

T' = ifK[fBn( m K)( m B - m l) + IbA^kWr}^ 

A' = ifB{fUO{™ 2 K - K) + fK,(m 2 B )m 2 B ]. (35) 

The amplitudes U and V have the same form as U and V, respectively, except that the 
CP-conjugate amplitudes contain the complex conjugate CKM matrix elements V u bV* s and 
V t bV t * s . As expected the strong and dipole penguin operators do not contribute to U and V. 
The only difference between the amplitudes (U and V) and the CP-conjugate cmplitudes 
(U and V) is the opposite sign of weak phase angle 7. We note that if electroweak penguin 
contributions are neglected, we would find that the relation of eq.(p2|) holds. However, we 
now obtain the following ratios 

\U - Ue 2i ^\ 



\U< 



166%|sin 7 | (36) 



tree 



and 

\V - Ve^l 



42% I sin 7 |, (37) 



I ^tree \ 

where the amplitudes U tree and V tree denote the tree contribution of U and V , respectively 
Here we have used m s = 175 MeV, fx = 162 MeV, fs = 200 MeV, and the form factors 



calculated in Ref. [p3| . We conclude that the assumption presented in Ref. is invalid and 



so the suggested isospin analysis is unworkable. 

IV. MEASUREMENT OF THE ANGLE (3 

In this section we study penguin effects on the determination of the CKM phase (3. Many 
methods have been suggested involving B meson decay into charmed particles. We consider 
two of the most convenient experimentally. 

15 



(1) Using B -> ^>K S 

The easiest way to measure j3 is to measure the parameter ImX in CP asymmetry of 
time evolution in B° —>■ ipKs 0. In this case, 

pA(B° ^ipK s ) J K ' 

The tree amplitude is proportional to V cb V* s , and the dominant penguin contribution 
from the internal top quark is proportional to V tb V t * s , and so the decay amplitude can be 
parametrized as 

A(B° - $K S ) = VaVZT+K + V tb V*P^ K . (39) 

Using unitarity of the CKM matrix, we can rewrite the above as 

A(B° -> $K S ) = V cb V*(T fK - P^k) - V ub V: s P^ K . (40) 

The WC's involved indicate that \T^ K \ is much larger than \P^k\- Also |KbKsl ^ s aDou t 
50 times larger than |V^6V^ S | from experimental data. We can, then, safely negelect the 
contribution from the term proportional to V ub V* s . To a very good approximation even if the 
penguin (strong and electroweak) effects are included, ImX = Im((q/p)(V c hV* s /V* b V cs )) = 
— sin(2/5). (3 can be measured with an error less than a percent. 

(2) Using B -> D+D- decay 

The same conclusion can not be drawn for method to determine (3 by measuring the 
parameter ImX in the process B° — ► D + D~ . In this case, the decay amplitude can be 
parametrized as 

A(B° -> D+D') = V cb V c * d T DD + V tb V t * d P DD 

= V cb V; d (T DD - P DD ) - V ub V* d P DD . (41) 

In this case although the penguin amplitude Pdd resulting from the top loop is suppressed 
compared with the tree amplitude Tdd, the CKM elements involved are comparable for 
each term. The error caused by penguin effects are much larger. If we keep the pegnuin 
contribution, we find 
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j x = sin(2/3) + 2R DD sin(3/3) cos 5 + R 2 DD sin(4/3) 

l+i% + 2IiWos(2/3 + <5) ' 1 j 

where Rdd = \Pdd/Tdd\, and 5 is the relative strong rescattering phase between the tree 
and penguin amplitudes. 

In factorization approximation we can calculate the ratio Rdd- For the decay B° — > 
D+D-, we find 

- D+u-) = ^{y c6 ^(eci + cso 

- KjVSKca + c 4 + 2(£c 5 + c 6 + £c 7 + c 8 )X' + (£c 9 + c 10 )]T', (43) 



where 



x>- m ° 



(m b - m c )(m c + m d ) ' 
= if D [f£ D (m 2 D )(m 2 B - ml) + /bd«)^]- (44) 

Using the effective coefficients q, the masses m D = 1.869 GeV, m c = 1.3 GeV, the decay 



constant fjj = 162 MeV, and the form factors calculated in Ref. |23|, we obtain 



Rdd = 0.09 . (45) 

Here we have neglected a small contribution due to a u-quark loop in penguin diagram 
(about 8% compared with a top-quark loop contribution) and annihilation effects. In Fig.3. 
we plot Asin(2/3) = 7mA + sin(2/3) as a function of (3. For the strong rescattering angle 5, 
we use the quark level estimate 5 = 12.4° by including absorptive contribution in the WC's. 
The error for certain values of j3 can be quite large. For example, for (3 = 45°, the error is 
above 16%. We also carried out calculations with the dipole penguin operator contribution. 
We again find their effects to be small. 



V. MEASUREMENT OF THE PHASE 7 



In this scetion we comment on some methods for measuring the CKM phase angle 7. 
Several different classes of method to measure 7 have been proposed. One class involve 
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charged B decays into a charge kaon and a neutral charmed meson [p7j . This class of 



methods is not contaminated by penguin effects because only tree amplitude can mediate 
such decays. These should work well and we do not consider them here. Another class of 
method is to use information from B — > mi, and B — > ir(r])K. In this case both tree and 
penguin amplitudes contribute, and care must be taken to include pegnuin effects. 

A method to measure 7 using these decays was first studied by Gronau, London and 
Rosner (GLR) [fTo] . In this method the tree and strong penguin effects were considered, but 
the electroweak penguin effects were neglected. GLR argued from isospin analysis that the 
combined amplitude A(B~ — > tt~K°) + A(B~ — > tt°K~) is a pure / = 3/2 amplitude, and 
therefore only the tree amplitude contributes. Using SU(3) relation shown in eq.fllOD, this 



amplitude is found to be equal to the decay amplitude (/i<-// 7r )(V u * s /\4rf) 2 A(i? _ 



7T 7T 



which can be measured. These three amplitudes form a closed triangle, and similarly for 
their CP-conjugate amplitudes. Using the fact that \A(B~ -> vr"^ )! = \A(B + -> 7r + fr°)| 
(because the tree contributions here are negnigiblely small), from the two triangles for the 
particle and anti-particle decay amplitudes, the relative phase 27 between A(B~ — > tt~7t°) 
and A(B + — > 7r + 7r°) can be obtained. This is a very interesting proposal. However it 
was soon pointed out by Deshpande and He [T^] that the inclusion of electroweak penguin 
effects invalidate this method because the electroweak contributions to / = 3/2 amplitude 
are compareable to the tree contribution. 

Other methods have been proposed to take into account the electroweak penguin effects. 
Recently Gronau, Hernandez, London and Rosner (GHLR) [13| showed that the difficulty 
with the electroweak penguin effects can be solved by constructing the quadrilaterals dis- 
cussed in section III for B — > ttK decays. As is already shown, the way used to construct the 
quadrilaterals in Ref. |J is not workable. Instead, GHLR used SU(3) relation to relate one 
of the diagonal of the quadrilateral to Bj — ► K~r\. This time the common side of the two 
quadrilaterals is chosen to be |v4(-B + — > n + K )\ = \A(B~ — ► 7r~^T°)|. This method is how- 
ever very difficult to implement experimentally because the decay amplitude for B~ — > K~rj 



is dominated by electroweak penguin contribution and has a very small branching ratio |28 



A more practical method has recently been proposed by Deshpande and He Jl4| using 
SU(3) relations between the decay amplitudes for AS = 1 decays B~ — > tt~K°, n°K~ , r]K~, 
and AS* = decay B~ — > 7t~tt° . This method requires the construction of the triangles 
obtained from eq. (ITT]) 



V2A{B~ -> ir°K-) - 2A{B~ -> -r-R ) = VgA(B~ -> r^fsT), 

V^A(S + -> 7T°^ + ) - 2A(5+ -> 7T+K ) = \/6A( J B+ -> 7? 8 if + ), (46) 

where r/g is the pure octet component. Using the relation 

A(B- -> tt-K°) = A(5+ -> 7T+K ), (47) 



the following result was obtained in Ref. [14] 



B-B = -t2V2e i5T ^4 1 7r_7r °) I sin 7' ( 48 ) 

I Vw| 

where £> and 5 are the complex quantities defined as 

B = V2A{B- -> vr ^-) - -> 7T-K ), 

B = V2A(B + -> 7T°^ + ) - A(5+ -> 7T + i?°), (49) 

shown in Fig.4. The angle 5 T denotes the strong final state rescattering phase of the tree 
amplitude of B (or B). Thus sin7 can be determined from eq. fl48|) . This method is free 
from the electroweak penguin contamination problem, and all decays involved have relatively 
large (O(10~ 5 )) branching ratios. They are within the reach of future experiments . 

The results given in Ref. [T4J hold in the exact SU(3) limit. This relation may be broken 
by SU(3) breaking effects due to 77-77' mixing, the breaking effects in form factors and mass 
differences. We now make quantitative estimates of the influence of SU(3) breaking effects 
in the factorization approximation. In this approximation we find 

Fx = V2A(B- -> tt°K-) - 2A(B- -> tTK ) 
= ^{V ub V: s [( Cl + £c 2 )C" + (fo + c 2 )(T' - A')] 
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- V tb V*[-^c 3 + c 4 )(T' + A') - 2(£c 5 + C6 )(YT' + ZA') 
+ 2(£c 7 + c 8 )(2FT' - ZA') -\{c 7 + £c 8 -eg- £c 10 )C" 
+ ^c 9 + c 10 )(2T' -A')}}, 

F 2 = V6A(B- -> 

= ^{K.C[(ci + £c 2 )C" + (fo + c 2 )(T" - A")] 

- V tb V*Mc 3 + c 4 )(T" - 2C" - A") + 2(£c 5 + c 6 )(-2X"C" + YT" - ZA") 
+ 2(£c 7 + c 8 )(X"C" + YT" - ZA") - ^(c 7 + £c 8 - c 9 - £c 10 )C" 

+ (Zc 9 + c 10 )(C" + T"-A")}}, (50) 

where C, T', A', Y and Z are given in eq. (p5|) and 



2m s (m b - m s ) ' 

C" = if m [fBK( m2 r, 8 )( m B ~ m \) + fB K ( m2 m ) m vs\> 
T" = ifK[fB Vs ( m K)( m B ~ m l s ) + fB m ( m K) m K], 

A" = if B [ti m (m 2 B)( m K - ml) + IkJOO- (51) 

We note that in SU(3) limit the triangle relation eq. (Hp) is verified. In numerical estimates, 
neglecting small contribution of the annihilation diagram, we obtain 

F x = ^i[V ub V: s (1.86GeV 3 )e i5T + V tb V*(-4.82 x 10- 2 GeV 3 )e l5p ] , 

F 2 = -}=i[V ub V* s (2.7QGeV 3 )e i5T + V^* s (-4.79 x lQ- 2 GeV 3 )e i6p ] . (52) 
v2 

In the above we have inserted arbitrary strong rescattering phase in the amplitudes. We 
can estimate the phases using absorptive part in WC's which indicate small phase for Sp 
and zero phase for St- We however keep them as free parameters here for covenience. For 
our numerical values we have used the decay constants f m = 176 MeV, f Vo ~ f m ( rj is 



the singlet component), and the from factors obtained in Ref. [ 2"3fl . We see that the SU(3) 
breaking effects in the tree amplitude are about 30%, and much smaller effects in the penguin 
amplitudes. 
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To have an idea how large the SU(3) breaking effects on the determination of 7 are, we 
carried out an excercise by taking i*\ 2 in eq.(|5^) to be the experimental values keeping 7 
and 5 = 5p — 5t as free parameters. For given values of 7 and 5, we construct a triangle and 
obtain the value for \B — B\. We then take the calculated amplitude for \A(B~ — > 7r~7r°)| 
as the measured value, and use eq.(fI8D to determine | sin 7). For each given 7, using eq.(fI8D 
we will obtain a output 7. We will call it 7'. Corresponding to the cases 7 > 5 and 7 < 5, 
there are two solutions for B — B which arise from the two possible orientations of the two 
triangles relative to their common side. Fig. 4. shows the two triangles used to find the 
magnitude \B — B\ for 7 > 5. We expect that because of the SU(3) breaking effects the 
triangle relation will have some deviation. This deviation will cause errors in determining 
I sin 7 1 and the angle 7. In Fig. 5. we show the errors A7 = 7 — 7' for a fixed strong phase 
5 = 12° which we expect from quark level evaluation of absorptive parts. Since 6 is expected 
to be small we focus on the case 7 > 5. From Fig. 5 we see that there are limited range 
where there is solution for 7' because a triangle can not be formed after breaking effects for 
all given 7 and 5. We see that the errors A7 increases as 7 increases. For instance, A7/7 is 
within about 20% for 23° < 7 < 35°. For larger values of 7 the error is larger. If the form 
factors are varied by taking a different fit, it is possible to reduce errors. We conclude that 
this method is very sensitive to SU(3) breaking effects. More theoretical efforts to study 
SU(3) breaking effects are called for. We hope Lattice calculation will provide us with useful 
information on the evaluation of amplitudes. 



VI. CONCLUSIONS 

We have studied the influence of penguin (especially, electroweak penguin) contributions 
for several methods to extract the angles a, j3, and 7 of the CKM unitary triangle. Our 
calculations are based on the factorization approximation using the general effective Hamil- 
tonian to the next-to-leading order and on some models for form factors when we have 
needed to obtain numerical values. To see the sensitivity of the results on the form factors 
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used, we repeated our calculations using form factors obtained in Ref. p0| . The conclusions 
are not changed for all cases except the calculation in Sec.V. Here the SU(3) breaking effects 
become larger. 

In the cases of measuring the phase a using isospin relations for B — ► tttt [|7| and sim- 
ilarly for the method using time dependent asymmetries in B — > pn ||, even though the 
electroweak penguin contributions contaminate the results, the phase a can be determined 
with an accuracy better than a few percent because the electroweak penguin effects are 
found to be quite small. However, for the approach proposed in Ref. |Tl|, we have shown 
that the method is unworkable since in this case the electroweak penguin effects are com- 
parable to the strong penguin effects. Similarly for method proposed in Ref. || to extract 
a from analysis of B — > ttK, we found that the electroweak penguin contributions can not 
be neglected and the assumption of the analysis is again invalid. 

For the measurement of /3, the penguin effects are negligible when use B — > ipKs- We 
found that penguin effects are not so small in extracting sin (2/3) from measurement of a 
CP-asymmetry in B — > D + D~. The deviation from sin (2/3) due to penguin contributions 
will be over 10% if the value of (3 is in the range of 12° < (3 < 62°. 

Finally we have made quantitative estimates of SU(3) breaking effects for the method of 
measuring the angle 7 proposed in Ref. Ijj . We found that the results are very sensitive to 
SU(3) breaking effects and permit extraction of 7 only in a limited range of parameter space. 
More study of the SU(3) breaking effects on the decay amplitudes are called for. Recently, 



methods using nonet symmetry [jn| have been suggested. These methods are subject to 
nonet symmetry breaking, and detailed studies are required to test their feasibility. 
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FIGURES 
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FIG. 2. The two quadrilaterals with a common diagonal U + V. Lines a, b, c, and d denote the 
amplitudes A(B° -» vr ^ ), A(B° -> 7r-^T + ), 4(5+ -» vr ^), and -> vr+iT )- Similarly 

lines a', b', c', and d! denote the corresponding amplitudes A's. 
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FIG. 3. A sin (2(3) = ImX + sin (2/3) versus (3 for a fixed strong rescattering phase 5 = 12.4°. 




FIG. 4. The triangles used to find \B — B\ for 7 > 5. Lines a, b, and c denote the amplitudes 
2 1 / 2 A(B- vr ^-), 1(5- -» tt-K°), and 6 1 / 2 A( J B~ -> r/gi^"). The dashed lines are for the 
corresponding B + decay amplitudes. 
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